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INTRODUCTION 


I . 


Carbon  materials  are  widely  used  for  thermal  protection  system 
applications  involving  very  high  temperatures  and  energy  fluxes.  The 
objective  of  our  four-year  concracc  research  program  has  been  to  gain  a  better 
understanding  of  the  superior  energy-dissipating  efficiency  of  these  materials 
and,  in  particular,  to  determine  the  limits  of  their  outstanding  performance. 
During  the  first  two  years  of  effort,  three  uncertain  aspects  of  carbon 
ablation  behavior  were  studied: 

a.  Vapor  pressure  of  carbon 

b.  Melt  temperature  of  carbon 

c.  Vaporization  kinetics  of  carbon 

Carbon  laser  vaporization  experiments  and  a  nonequil ibrium  ablation 

analysis  were  carried  out  to  obtain  basic  information  concerning  the  above 

properties.  Our  first  interim  report^  contains  a  description  of  the 

experiments,  the  analytical  approach,  and  a  preliminary  analysis.  Detailed 

data  analyses  and  conclusions  concerning  these  carbon  thermochemical 

2 

properties  are  given  in  the  second  interim  report,  as  well  as  in  Refs.  3 
and  4  . 


^Baker,  R.  L. ,  "Carbon  Nonequilibrium  Phase  Change,"  Office  of  Naval  Research 
Interim  Report,  Tr'- 0081(6728-02 )-l ,  The  Aerospace  Corporation,  El  Segundo , 

CA,  Dec.  1981. 

^Baker,  R.  L.  and  M.  A.  Covington,  "The  High  Temperature  Thermochemical 
Properties  of  Carbon,"  Office  of  Naval  Research  Interim  Report, 

TR-0082 ( 2729 )-l ,  The  Aerospace  Corporation,  El  Segundo,  CA,  Mar.  1982. 

^Baxer,  R.  L. ,  M.  A.  Covington,  and  G.  M.  Rosenblatt,  "The  Determination  of 
Carbon  Thermochemica I  Properties  by  Laser  Vaporization,"  High  Temperature 
Materials  Chemistry  Symposium,  The  Electrochemical  Society  (1983)  pp.  143-134. 

^Baker,  R.  L. ,  M.  A.  Covington,  and  G.  M.  Rosenblatt,  "The  Vapor  Pressure, 

Melt  Temperature  and  Vaporization  Kinetics  of  Carbon,"  J .  App 1 .  Phy s . , 

(To  be  published.) 
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During  our  third  year  of  effort,  laser  ablation  experiments  and  analyses 
were  carried  out  and  reported  for  metal-impregnated  carbon  materials.  In 
addition,  higher  order  numerical  simulations  of  a  nonequilibrium  ablation 
process  were  performed  and  reported. This  latter  work  confirmed  the 
approximate  analysis  methods  we  had  previously  used  for  obtaining  carbon 
thermochemical  properties  from  our  laser  ablation  data. 

This  final  report  describes  additional  experiments  and  analyses  carried 
out  during  the  past  year.  Our  present  objectives  are  closely  related  to  prior 
work.  They  include  determination  of  both  the  relative  surface  emissivities 
and  the  vaporization  kinetics  of  different  types  of  graphite.  Dnissivity  and 
vaporization  kinetics  effects  are  discussed  in  Sections  II  and  III, 
respectively. 

An  additional  objective  of  recent  work  has  been  to  investigate  other 
effects.  It  is  of  particular  interest  and  importance  to  determine  the 
character  and  extent  of  carbon  particulates  in  the  expanding  carbon  vapor  and 
their  effect  on  the  transport  of  infrared  wavelength  radiation  through  the 
plume.  A  new  test  chamber  and  infrared  spectroscopic  instrumentation  were 
fabricated  tnis  past  year  to  study  these  effects.  The  test  apparatus, 
instrumentation,  and  initial  experimental  results  are  discussed  in  Section  IV. 

An  overall  summary  and  conclusions  from  our  four-year  effort  are  given 
in  Section  V. 


^Turner,  J.  S. ,  "Laser- Induced  Phase  Changes  at  Solid  Surfaces,"  Final  Report, 
The  University  of  Texas  at  Austin,  October  1982. 

^Baker,  R.  L. ,  D.A.  Nelson,  and  J.  S.  Turner,  "Higher  Order  Numerical 
Simulations  of  the  Knudsen  Layer,"  Office  of  Naval  Research  Interim  Report, 
TR-0083(3729)-2,  The  Aerospace  Corporation,  El  Segundo,  CA,  Mar.  1983. 
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II.  EMISSIVITY  DATA  AND  ANALYSIS 


Our  interest  in  determining  the  high  temperature  emissivities  of 
different  graphites  stems  from  the  vaporization  kinetics  studies  we  have 
carried  out.  As  shown  in  Refs.  2-4,  and  further  discussed  in  Section  III,  the 
experimentally  measured  relative  vaporization  kinetics,  for  different 
graphites,  are  very  sensitive  to  the  values  of  the  surface  emissivity  used  to 
determine  the  surface  temperatures.  Because  of  this  and  a  lack  of  high 
temperature  emissivity  data  for  ablation-roughened  samples  in  the  literature, 
new  experiments  were  specifically  designed  and  carried  out  to  determine  this 
information.  The  experiments,  data,  and  analysis  are  discussed  below. 

A.  EXPERIMENTAL  TECHNIQUE 

1 .  Approach 

Due  to  difficulties  associated  with  making  measurements  ac  temperatures 
above  2000°K.  and  the  small  ablated  (roughened)  sample  areas  available,  it  was 
necessary  to  design  an  experimental  technique  specifically  for  our  purposes. 
Dr.  C.M.  Rosenblatt*  suggested  that  small  samples  of  the  graphite  for  which 
Che  emissivity  is  desired  could  be  embedded  in  the  surface  of  a  bulk  sample  of 
a  graphite  for  which  the  emissivity  is  known.  Then,  by  bringing  the  bulk 
graphite,  with  its  embedded  sample,  up  to  a  high  temperature  and  measuring  the 
irradiance  from  its  surface  and  from  the  sample  surface,  the  emissivity  of  the 
cample  could  be  determined  from  the  relative  irradiance  values.  If  it  is 
assumed  chat  the  temperatures  of  the  two  materials  are  the  same,  then  the 
relative  irradiance  is  a  direct  measure  of  Che  relative  emissivities. 


*  Personal  communication,  Los  Alamos  National  Laboratory,  Los  Alamos,  N.  M. 
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2. 


Apparatus 


As  originally  envisioned,  the  bulk  graphite  was  to  be  resistively 
heated.  For  the  experiments  that  were  carried  out,  the  graphites  were 
laser-heated.  The  "bulk"  graphite  sample-holder  consisted  of  a  pyrolytic 
graphite  cylinder  8.5  mm  in  diameter  and  4  mm  long  with  a  3-mm  diameter  sample 
hole  on  its  axis.  A  photograph  of  two  saraple-holder9  is  shown  in  Fig.  la.  A 
second  hole,  shown  in  Che  cylindrical  surface  of  the  holder  on  the  right, 
accommodated  a  post  used  to  hold  the  sample-holder  and  sample  in  place.  A 
photograph  of  the  entire  assembly,  including  sample-holder,  sample,  and  post, 
positioned  inside  a  5-mra  thick  graphite  radiation  shield  is  shown  in  Fig.  lb. 

In  order  to  irradiate  one  face  of  the  sample  and  holder  with  a  laser 
and  measure  the  relative  irradiauce  on  the  opposite  face,  the  radiation  shield 
with  the  sample  assembly  inside  was  attached  to  a  1/2-in.  thick  graphite  plate 
and  positioned  onto  the  centerline  of  the  vacuum  chamber  apparatus  described 

in  Section  V.  This  setup  is  shown  schematically  in  Fig.  2.  The  laser  beam 

was  defocussed  in  order  to  uniformly  flood-load  the  face  of  the  sample  and 

holder.  Using  a  quartz  viewing  window  located  in  a  port  on  the  rear  of  the 

chamber,  we  found  it  possible  to  directly  view  the  opposite  face  with  a 
scannable  pyrometer. 

A  close-up  phocograph  of  a  sample  and  holder  in  the  radiation  shield  is 
shown  in  Fig.  3a.  The  pyrolytic  graphite  holders  were  fabricated  so  that  the 
planes  of  this  highly  anisotropic  material  were  perpendicular  to  the  axis  of 
the  support  post.  This  accomplished  several  important  objectives: 

•  It  minimized  heat  transfer  losses  into  the  supporting  post,  which  was 
also  made  from  pyrolytic  graphite  with  the  planes  perpendicular  to  the 
post  axis.  (The  thermal  conductivity  perpendicular  to  the  planes  is 
200  times  less  than  in  the  directions  parallel  to  the  planes. 


^Touloukian,  Y.  S.,  ed.,  Thermophysical  Properites  of  Matter,  Vol.  8, 
IFI/Plenum  Press,  NY,  1970. 
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Fig.  1.  Emiss 


Fig.  2.  Sample  and  Holder  Mounted  for 
Laser  Irradiation 


'  It  insured  a  very  high  thermaL  conductivity  from  the  front  (laser 
irradiated)  to  the  back  of  the  sample-holder,  thus  minimizing 
temperature  gradients  and  maximizing  laser  energy  transport  into  the 
sample-holder. 

*  It  exposed  the  edges  of  pyrolytic  graphite  planes  to  both  the  incoming 

laser  irradiation  as  well  as  the  viewing  pyrometer.  Since  the 
7  S 

literature  ’  supports  a  unit  emissivity  value  for  pyrolytic  graphite 
looking  edge-on  to  the  planes,  the  absorbed  laser  energy  was  maximized 
and  the  relacive  emissivity  of  the  sample  in  the  center  could  be 
compared  with  the  unit  emissivity  of  the  sample-holder. 

The  surfaces  of  the  3-mm  diameter  x  4-mm  long  cylindrical  samples,  to 
be  viewed  by  the  pyrometer,  were  pre-roughened  by  focussing  the  laser  on  them 
for  1-2  sec  in  an  atmospheric  air  environment.  A  photograph  of  roughened 
samples  of  the  three  graphites  for  which  emissivity  measurements  were  made  is 
shown  in  Fig.  3b. 

3 .  Operating  Procedure  and  Data 

The  1  kW  CO^  laser  used  to  heat  the  sample-holder  with  a  sample 
inserted  is  described  in  Section  V.  Normal  operating  procedure  called  for 
irradiation  at  nominal  on-sample  laser  power  levels  of  500,  700,  and  900  W. 

At  each  power  level,  sufficient  time  (15-30  sec)  was  allowed  for  steady-state 
temperatures  to  be  reached.  The  pyrometer  was  then  mechanically  scanned 
horizontally  across  the  entire  sample-holder  along  a  line  passing  through  the 
centerline  (axis)  of  the  sample.  The  pyrometer  (Thermogage  Model  8000-10)^ 
had  a  nominal  spot  size  of  1/16  in.  (1.59  mm)  and  a  focal  length  of  30  in. 

(76  mm).  Remotely  operated  and  electrically  driven  micrometers  were  used  to 


8£z  rolytic  Graphite  Handbook,  Metallurgical  Production  Department,  General 
Electric  Company,  Schenectady,  NY,  1964,  as  quoted  in  Gokcen,  N.A, ,  et.  al., 
"Determination  of  Graphite/Liquid/Vapor  Triple  Point  by  Laser  Heating,"  High 
Temperature  Science,  Vol.  8,  June  1976,  pp  31-97 

^Brookley,  C.E.,  "Product  Data  for  High  Temperature  Pyrometers,"  Thermogage, 
Inc.,  Krostburg,  MD,  1982.  Also,  "Calibration  Curves  for  the  Aerospace  High 
Temperature  Pyrometer." 


sweep  the  pyrometer  spot  from  one  edge  to  the  other  edge  of  the  sample-holder 
in  about  30  sec.  After  two  or  more  scans  at  a  given  laser  power  level,  the 
power  was  increased  and  the  procedure  repeated. 

Experimental  data  obtained  in  this  manner  for  graphite  samples 
consisting  of  Graphnol  N3M^ ,  Poco  (AXF.Q1),  and  pyrolytic  graphite  (Super 
Temp,  continuously  nucleated)  are  shown  in  Figs.  4,  3,  and  6,  respectively. 

In  each  figure,  the  relative  sizes  of  the  sample-holder,  sample,  and  pyrometer 
spot  are  shown  at  the  top;  the  pyrometer  signals  obtained  by  scanning  at  each 
of  chree  laser  power  levels  are  seen  in  the  lower  half  of  the  figure. 

B.  DATA  ANALYSIS 

9 

When  c  inverted  to  temperatures,  the  pyrometer  output  signals  (mV), 
shown  in  Figs.  4,  5,  and  6,  give  values  ranging  from  1780°K  to  2120°K  for  a 
surface  emissivity  of  1.00.  For  emissivity  e  less  than  1.00,  the  tempera¬ 
ture  is  obtained  from  the  requirement  that  the  surface  irradiance  is  a 

_4  „ 

constant;  i.e.,  £T  =  constant,  for  a  given  data  point. 

The  pyrometer  data  shown  in  Figs.  4  and  5,  for  the  polycrystalline 
graphites  Graphnol  and  Poco,  respectively,  indicate  a  relatively  flat  output 
voltage  profile  across  the  samples.  As  the  laser  power  is  increased,  higher 
voltages  (temperatures)  are  obtained;  but  the  profiles  across  the  samples 
remain  reasonably  flat.  Based  on  our  assumption  in  designing  the  experiment, 
that  the  temperature  of  the  sample  is  the  same  as  that  of  the  sample-holder, 
we  conclude  from  these  figures  that  the  roughened  surface  emissivities  of 
Graphnol  and  Poco  graphites  are  nearly  the  same  as  the  unit  emissivity  of  the 
flat  surface  of  the  pyrolytic  graphite  sample-holder. 

l°Johnsen,  B.  P.  and  H.  S.  Starrett,  "Preliminary  Design  Data  Package  for 
Graphnol  N3M  (Phase  III,  Mantech),"  Report  SORI-MER-82-20 ,  Southern  Research 
Institute,  Birmingham,  AL,  Nov.  1982. 
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In  contrast  to  che  above,  the  data  shown  in  Fig.  6  for  the  pyrolytic 
graphite  sample  indicate  substantially  less  output  signal  from  the  pyrometer 
when  it  is  scanning  the  sample.  If  the  sample  emissivity  is  1.00,  then  these 
data  indicace  sample  temperatures  110  to  12Q°K  lower  than  the  temperatures  of 
the  sample-holder.  However,  a  more  likely  explanation  is  that  the  sample  and 
sample-holder  temperatures  are  nearly  the  same,  for  a  given  laser  power  input 
level,  and  the  lower  pyrometer  output  signal  is  due  to  a  reduced  emissivity  of 
the  sample  relative  to  the  sample-holder.  Interpreted  in  this  way,  the  data 
in  Fig.  6  indicate  a  surface  emissivity  of  0.80  for  the  pyrolytic  graphite 
sample  viewed  perpendicular  to  che  lamellar  planes.  This  should  be  regarded 
as  a  lower  limit  for  the  emissivity  of  this  material  in  this  direction  as 
discussed  below. 

As  previously  mentioned,  pyrolytic  graphite  has  a  thermal  conductivity 
parallel  to  the  lamellar  planes  which  is  200  times  greater  than  that 
perpendicular  to  the  planes.  For  the  Fig.  6  data  discussed  above,  laser 
energy  was  conducted  through  the  pyrolytic  graphite  sample-holder  along  the 
planes.  In  contrast,  conduction  through  the  pyrolytic  graphite  sample  was 
perpendicular  to  the  planes.  This  leads  to  obvious  questions  concerning  our 
assumption  that  the  sample  and  sample-holder,  as  viewed  by  the  pyrometer  (see 
Fig.  2),  are  at  the  same  temperature. 

Less  than  one  millimeter  below  the  surface  where  the  laser  energy 
enters,  the  sample  temperature  is  the  same  as  that  of  the  surrounding  sample- 
holder.  This  conclusion  is  based  upon  a  highly  simplified  two-dimensional 
model  of  che  energy  transport  and  temperature  distributions  through  this  sample 
and  sample-holder.  From  this  point  until  the  back  surface  is  approached,  very 
little  energy  is  transported  through  the  9ample  along  its  axis  because  of  the 
extremely  low  thermal  conductivity  in  this  direction.  At  each  axial  location 
beyond  this  paint,  the  sample  temperature  is  the  same  as  that  of  the  surround¬ 
ing  sample-holder  due  to  very  high  sample  thermal  conductivity  in  the  radial 
direction.  As  the  pyrometer-viewed  (back)  surface  is  approached,  the  sample 
temperature  drops  rapidly  in  order  to  conduct  energy  to  the  surface  to  be 
radiated  away.  The  approximate  energy  balance  equations  in  this  region 
indicace  that  the  temperature  seen  by  the  pyrometer  on  the  sample  centerline 
could  be  30°K  less  than  che  temperature  seen  on  the  sample-holder. 
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If  this  estimate  is  reasonable,  the  drop  in  pyrometer  output  signal 
(seen  in  Fig.  6  when  the  sample  is  viewed)  is  due  to  a  decreased  emissivity  of 
the  sample  relative  to  the  holder  and  to  a  sample  temperature  slightly  less 
than  Chat  of  the  sample-holder.  In  this  case,  rather  than  the  lower  bound 
emissivity  of  0.80  previously  found,  this  interpretation  of  che  data  indicates 
an  emissivity  of  approximately  0.85. 

An  alternate  method  of  determining  the  high  temperature  emissivity  of 
pyrolytic  graphite  has  also  been  developed.  Unfortunately,  within  che 
constraints  of  our  laser  operating  time,  we  were  not  able  to  obtain  additional 
data  to  test  this  new  method.  The  technique  and  approacn  are  discussed  in 
Appendix  A. 
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III.  GRAPHITE  VAPORIZATION  KINETICS 


An  important  major  conclusion  of  our  earlier  work  was  that  the 
vapor  pressure  of  carbon  is  "best  represented"  by  thermodynamic  property  data 

other  than  that  given  in  Che  JANAF  Tables.^  The  recommended  data  are  those 

12  13 

of  Lee  and  Sanborn  and  Leider,  Knkonan,  and  Young.  Since  these 

workers  were  at  the  Lawrence  Livermore  National  Laboratory  when  their  results 

were  published,  we  have  referred  to  these  data  as  the  Livermore  property 

data.  When  these  data  are  used,  the  predicted  carbon  total  vapor  pressures 

are  cwo  to  four  times  greater  than  those  obtained  using  Che  JANAF  property 

data. 

As  discussed  below,  the  rate  of  mass  loss  from  a  solid  carbon 
(graphite)  surface  is  proportional  to  both  the  vapor  pressure  and  the 
vaporization  coefficient  a.  Whereas  the  vapor  pressure  is  an  equilibrium 
property  of  Che  material,  the  vaporization  coefficient  is  a  nonequilibrium 
(kinetic)  parameter  which  expresses  the  relative  rate  of  (in  this  case)  the 
carbon  solid-vapor  sublimation  process.  Thus,  the  importance  of  obtaining 
carbon  vaporization  coefficient  conclusions  consistent  with  our  Livermore 
thermodynamic  property  data  conclusion  becomes  apparent. 

Preliminary  conclusions  regarding  numerical  values  of  carbon 
vaporization  coefficients  were  given  in  Refs.  2-4.  Insofar  as  we  know,  these 
are  the  only  publications  in  the  literature  that  give  an  experimentally 
determined  temperature  dependence  of  these  coefficients.  Subsection  A,  below, 
review'  basic  definitions  and  the  importance  of  vaporization  coefficients.  In 


1 ‘JANAF  Thermochemical  Tables,  National  Bureau  of  Standards,  NBS-37,  June 
19  71. 

l^Lee,  E.  L.  and  R.  H.  Sanborn,  "Extended  and  Improved  Thermal  Functions  for 
the  Caseous  Carbon  Species  C  j-C  j  from  298  to  10,000°K,"  High  Temperature 
Sc ience ,  Vol.  5,  1973,  pp.  438-433. 

^Leider,  H.  R. ,  0.  H.  Krikorian,  and  D.  A.  Young,  "Thermodynamic  Properties 
of  Carbon  Up  to  the  Critical  Point,"  Carbon ,  Vol.  11,  1973,  pp.  555-363. 
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Subsection  B,  our  prior  results  and  conclusions  concerning  carbon  vaporization 
coefficients  are  presented.  Finally,  in  Subsection  C,  best  estimate  numerical 
values  of  a  are  given  for  pyrolytic  graphite  and  the  polycrystalline 
graphite,  Graphnol. 

A.  BASIC  DEFINITIONS  AND  IMPORTANCE  OF  VAPORIZATION  COEFFICIENTS 


Physically,  the  vaporization  coefficient  a  for  a  given  material 
represents  the  relative  vaporization  rate  of  that  material  at  a  given 
temperature.  It  is  defined  in  terms  of  the  maximum  possible  vaporization 
rate,  i.e. 
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The  vaporization  coefficient  defines  the  relative  rate  of  vaporization;  thus 
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^Rosenblatt,  G.  M.  ,  "The  RoLe  of  Defects 
Antimony,"  Surface  Defect  Properties  of 


in  Vaporization:  Arsenic  and 
Solids ,  Vol.  5,  1976,  pp. 36-64. 
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If  Che  Langmuir  vapor  pressure  can  be  determined  experimentally  as  a  function 

of  T  ,  then  by  knowing  p  (T  )  one  can  determine  the  temperature 
s  3  3  2  4 

dependence  of  cn  from  the  data.  1 

The  discussion  above  assumes  chat  the  vaporization  process  takes  place 
into  a  vacuum  with  no  vaporized  atoms  or  molecules  backscatcered  to  the 
vaporizing  surface.  In  Refs.  2  through  4,  an  analysis  method  based  upon  a 
model  of  a  nonequil ibriura  vaporization  process  s  developed  which  allows  the 
Langmuir  vapor  pressure  p.  to  be  determined  from  experimental  data  taken  in 
non-vacuum,  i.e.,  finite  back-pressure,  environments.  Once  the  temperature 
dependence  of  p^  has  been  measured,  a  as  a  function  of  Tg  can  be 
determined  in  Che  same  manner  as  discussed  above.  This  is  illustrated  is 
Subsections  B  and  C  below. 

Why  is  a  knowledge  of  vaporization  coefficient  values  important?  The 
primary  reason  is  that  for  materials  such  as  carbon  (graphite),  which  vaporize 
subject  to  kinetic  constraints,  the  rate  of  mass  loss  (ablation)  for  a  given 
energy  inpuc  to  the  surface  cannot  be  predicted  unless  both  the  vapor  pressure 
p^  and  the  vaporization  coefficient  a  are  known. 

Another  reason  we  need  to  determine  the  vaporization  coefficients  of 
different  graphite  materials  is  to  compare  their  vaporization  kinetics  to 
their  solid  state  structure  (crystallinity)  and  to  each  other.  Generally 
speaking,  for  a  given  rate  of  energy  input  to  a  graphite  surface,  the  material 
with  the  lowest  vaporization  coefficient  will  achieve  the  highest  surface 
temperature.  This  is  because  the  low  vaporization  coefficient  value  implies 
slow  kinetic  rates;  to  increase  the  vaporization  rate  to  that  required  for 
maintaining  an  energy  balance,  a  higher  temperature  is  needed.  Following 
through  on  this  idea,  a  graphite  material  with  vaporization  coefficients  less 
than  those  of  other  graphite  materials  can  be  melted,  i.e.,  it  reaches  the 
melt  temperature,  at  energy  flux  levels  lower  than  those  required  for  the 
material  with  higher  vaporization  coefficients.  Because  melting  degrades  the 
ablation  efficiency,  differences  such  as  these  are  important  to  understand  and 


B.  PRIOR  RESULTS  AND  CONCLUSIONS 


The  carbon  laser  ablation  experiments  which  were  carried  out 
specifically  to  evaluate  vaporization  coefficients  for  different  types  of 
graphice  are  described  in  Ref.  1.  Both  continuous  wave  (CW)  as  well  as  pulsed 
laser  experiments  were  performed.  The  mathematical  model  developed  to 
describe  the  carbon  nonequilibrium  ablation  process  is  described  in  Refs.  1 
and  2.  Data  analysis,  using  the  model  to  obtain  numerical  values  of 
vaporization  coefficients  from  the  measured  data,  is  contained  in  Refs.  2-4. 
Our  vaporization  coefficient  conclusions  from  these  analyses  are  summarized  in 
figs.  7  and  8. 

Figure  7  illustrates  our  results  for  pyrolytic  graphite.  The  lower 
solid  line  shows  the  measured  temperature  dependence  of  the  vaporization 
coefficient  as  obtained  from  the  CW  laser  data.  For  these  data,  both  the 
laser  illumination/  vaporization  and  the  pyrometer  viewing  took  place  in  a 
direction  parallel  to  the  pyrolytic  graphite  planes.  In  this  direction,  an 

g 

assumed  surface  emissivity  e  of  unity  is  supported  by  the  literature. 

This  alleviates  uncertainties  associated  with  measured  surface  temperatures 
and  lends  confidence  to  the  vaporization  coefficient  values  shown  in  Fig.  7 
for  the  temperature  range  from  3300  to  4000°K. 

The  upper  solid  line,  obtained  from  pulsed  laser  data,  has  a  slightly 
higher  slope  and  extends  measured  vaporization  coefficient  values  from  4000°K 
up  to  4500<>K.  We  note  three  observations  concerning  this  line.  First  of  all, 
the  line  is  for  an  arbitrarily  assumed  surface  emissivity  of  0.97.  Second, 
the  extrapolated  numerical  value  of  a  approaches  unity  at  the  Livermore 
triple-point  (melt)  temperature  of  4765°K.  Finally,  3t  4000°K  both  solid 
lines  indicate  approximately  the  same  values  of  a. 
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Measured  Vaporization  Coefficients 
for  Pyrolytic  Graphite 


The  laser  illuraination/vaporization  and  pyrometer  viewing  for  the 
pulsed  laser  data  (upper  solid  line)  were  in  a  direction  perpendicular  to  the 
pyrolytic  graphite  planes.  In  this  direction,  there  is  much  greater 
uncertainty  regarding  the  values  of  the  surface  emissivity  e.  The 
sensitivicy  of  the  data  interpretation  to  the  value  assumed  for  the  surface 
emissivity  is  illustrated  by  the  two  dashed  lines  in  Fig.  7  which  are  for  C 
values  of  0.9  and  1.0.  Thus,  the  importance  of  determining  the  correcc  value 
of  c  is  noted  (see  Section  II). 

Additional  experimental  data  were  obtained  in  both  the  CW  and  the 

pulsed  laser  experiments  to  evaluate  vaporization  coefficients  for  the 

polycryscal line  graphite,  Graphnol.  Unfortunately,  the  CW  laser  data  were 

ultimately  determined  to  be  uninterpretable  because  of  complications  related 

2 

to  the  necessary  presence  of  oxygen  in  the  chamber.  However,  numerically 
determined  vaporization  coefficients  for  both  pyrolytic  graphite 
(perpendicular  to  the  planes)  and  Graphnol,  obtained  from  the  pulsed  laser 
data,  are  shown  in  Fig.  8.  It  was  assumed  that  both  materials  had  a  surface 
emissivity  of  0.9.  For  this  assumption,  measured  values  of  a  for  Graphnol 
are  a  factor  of  5  lower  than  for  pyrolytic  graphite  at  a  given  temperature. 

The  overriding  importance  of  knowing  the  surface  emissivity  values  is  again 
apparent.  In  this  case,  little  can  be  said  about  the  relative  vaporization 
kinetics  of  the  two  graphites  until  their  associated  surface  eraissivities  are 
known. 

In  Section  II,  numerical  values  of  high  temperature  surface  emissivity 
were  determined  from  our  recent  experimental  data  for  ablation-roughened 
samples  of  Graphnol  and  pyrolytic  graphite  (perpendicular  to  the  planes).  In 
Subsection  III-C  below,  these  values  are  used  to  reinterpret  our  original 
pulsed  laser  data.  This  eliminates  the  arbitrary  emissivity  value  assumptions 
associated  with  the  above  discussion  of  Figs.  7  and  8  and  allows  "beat 
estimate"  values  of  a  to  be  determined  for  these  materials. 


C.  BEST  ESTIMATE  NUMERICAL  VALUES  OF  a  FOR  PYROLYTIC 
GRAPHITE  AND  GRAPHNOL 


In  the  following  discussion,  vaporization  coefficients  and  emissivities 
are  designated  all,  ell,  and  cti,  e_L  for  pyrolytic  graphite  parallel  and 
perpendicular  to  the  planes,  respectively.  When  the  emissivity  value  ( el)  of 
0.85  (determined  for  pyrolytic  graphite  in  Section  II)  is  used,  the  vaporiza¬ 
tion  coefficient  vaLues  (a-L)  interpreted  from  our  pulsed  laser  data  are 
snown  by  the  solid  line  on  the  left  in  Fig.  9.  Also  shown  by  the  line  on  the 
right  are  the  values  of  all  obtained  from  our  CW  laser  data.  We  note  that  at 
4000°K,  the  pyrolytic  graphite  vaporization  coefficient  aj_  for  vaporization 
taking  place  perpendicular  to  the  planes  is  about  a  factor  of  three  lower  than 
all  for  vaporization  occurring  parallel  to  the  planes.  A  plausible  explana¬ 
tion  of  this  behavior  is  that  for  vaporization  that  occurs  parallel  to  the 
planes,  i.e.,  from  an  a-b  face,  etch  pits  formed  in  the  vaporizing  surface 
create  more  "kimes"  or  active  sites^  than  when  the  vaporization  takes  place 
perpendicular  to  the  planes,  i.e.,  from  a  c-face.  The  greater  number  of  kink 
sites  in  the  former  case  leads  to  more  rapid  vaporization  rates  within  the 
context  of  the  classical  "terrace-ledge-kink"  model  of  solid  vaporization 
processes . ^  ^ 

A  comparison  of  high  temperature  vaporization  coefficients  for 
pyrolytic  graphite  and  Graphnol  is  shown  in  Fig.  10.  Again,  the  values  shown 
were  obtained  from  our  pulsed  laser  data  using  the  values  of  surface 
emissivity  determined  for  these  materials  in  Section  II.  Recall  that  the  same 
daca  were  interpreted  assuming  that  both  materials  had  an  emissivity  of  0.9  in 
Fig.  8.  Figure  10  indicates  that  within  the  limits  of  our  ability  to 
determine  emissivities,  mass  loss  rates,  and  temperature,  these  two  materials 
differ  very  little  in  their  kinetic  rates  of  vaporization. 


^Rosenblatt,  G.  M. ,  "Evaporation  from  Solids,"  Treatise  on  Solid  State 
Chemistry,  Vol.  6a,  N.B.  Hannay ,  ed. ,  Plenum  Press,  N.Y.,  1976,  p.  199. 

^Hirth,  J.  P.  and  G.  M.  Pound,  "Condensation  and  Evaporation,"  Progress  in 
Materials  Science,  Vol.  11,  MacMillan,  New  York,  1963. 
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PYROLYTIC  GRAPHITE 


Fig.  9.  Best  Estimate  Vaporization  Coefficients  for 
Pyrolytic  Graphite 
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IV.  PARTICLE  EFFECTS 

1 

Our  earlier  experiments'-  dramatically  showed  the  significant  carbon 
particulate  condensation,  and  consequent  laser-beam  blockage  effects,  which 
can  occur  when  carbon  is  laser-vaporized  into  an  inert  gas  (argon) 
atmosphere.  The  purpose  of  the  work  described  herein  was  to  begin  to  quantify 
these  effects  for  ambient  conditions,  ranging  from  atmospheric  pressure  to 
near-vacuum  conditions,  by  using  infrared  spectroscopic  methods  to  probe  the 
carbon  vapor-particle  plume  flowfield.  In  Subsection  A  below,  a  brief 
overview  is  given  of  the  spec croscopic  methods  we  have  used.  A  more  detailed 
description  of  the  experimental  apparatus,  spectroscopic  methods,  and  basic 
operating  procedure  io  given  in  Subsection  B.  Initial  experimental  results 
and  observations  are  then  discussed  in  Subsections  C  and  D,  respectively. 

A.  PARTICLE  EFFECTS  SPECTROSCOPY 

Infrared  spectroscopy  at  various  positions  in  a  C02  laser-generated 
carbon  plume  is  capable  of  providing  the  following  plume  diagnostic 
information : 

o  Spectral  transmission  profiles 
o  Temporal  variations  in  transmission 
o  Measurement  of  particulate  growth  rate 
o  Estimates  of  average  particulate  size 
o  Observation  of  molecular  constituents  (e.g.,  ) 

o  Laser  blockage  effects 

o  Plume-scattering  properties  at  10.6  pm 
o  Particulate  heating/cooling  rates 
o  Condensation  rates  in  vacuum  plumes 
o  Inert  gas  effects  on  condensation  rates 

Because  of  the  large  number  of  unknowns,  our  initial  efforts  can  be 
described  as  survey  experiments.  Low  resolution  infrared  spectroscopy,  using 
circular  variable  interference  filters  (CVIF),  was  selected  as  the  optimum 
method  of  performing  the  initial  diagnostics.  CVIFs  are  easy  to  use,  are 
capable  of  reasonably  fast  spectral  scan  time,  and  have  high  optical 
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through-put.  For  our  setup,  this  meant  1  to  2  percent  spectral  resolution  and 
a  0.75-sec  scan  time. 

Available  CVIFs  were  used,  which  restricted  our  measurements  to  a 
4-8  Um  scan  or  a  7-12  Um  scan.  A  typical  background  spectrum  is  shown  in 
Fig.  11.  Note  that  under  the  conditions  described  above  and  using  a  sensitive 
Hg:Cd:Te  detector,  cooled  to  77°K,  a  signal  to  noise  of  approximately  20/1  was 
available  for  these  studies.  This  provided  useful  spectroscopy  as  evidenced 
by  the  appearance  of  the  strong  6  )jn  H^O  absorption  from  the  external  1.2  m 
optical  path.  Certainly,  good  pLume  transmission  data  could  be  gathered  over 
large  portions  of  the  4-12  pm  spectral  region.  This  range  also  encompasses 
the  interesting  5.2  pm  rotation-vibration  band  for  the  molecule.^ 

B.  EXPERIMENTAL  APPARATUS  AND  OPERATING  PROCEDURE 

Our  experiments  were  conducted  in  the  2-ft  diameter  vacuum  facility 
shown  schematically  in  Fig.  12  and  pictorially  in  Figs.  13,  14,  and  15.  The 
major  facility  subsystems  comprise  the  vacuum  chamber,  CO^  laser  transfer 
optics,  and  the  diagnostics. 

1 .  Vacuum  Chamber 

A  stainless-steel  vacuum  chamber  was  available  and  measured  2  ft  in 
diameter  by  2.5  ft  in  length.  The  chamber  was  modified  to  focus  the  CC>2 
laser  energy  on  the  end  of  a  3-mm  diameter  x  25-mm  length  carbon  rod.  The 
rods  can  be  positioned  at  various  distances  from  the  front  of  the  chamber. 
Optical  ports  were  cut  in  the  sides  of  the  chamber  to  permit  infrared  spectral 
probing  through  KCJL  windows  across  the  effluent  carbon  plume.  The  ability 
to  reposition  the  carbon  rod  allowed  the  spectral  probe  diagnostics  to  remain 
stationary,  while  particle  effects  were  observed  at  various  distances  from  the 
laser-heated  carbon  source. 

^Treffera,  R.  R.  and  D.  P.  Gilra,  "The  Vibrational  Spectrum  of  C3  in  the 
5  Micron  Region,"  Astrophysical  J. ,  Voi.  202,  1975,  pp.  839-843. 
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Fig.  12.  Vacuum  Chamber,  Laser  Transfer  Optics,  and  Diagnostics 


Fig.  13.  Vacuum  Chamber  and  Optical  Bench 


■sls'i 


The  CO,  laser  beam  is  introduced  through  a  6-in.  diameter  port 
located  in  the  front  door.  The  6-in.  diameter  diagnostic  ports  are  mounted  in 
the  chamber  sides  on  a  horizontal  axis  perpendicular  to  the  laser  optical 
axis.  Additional  instrumentation  can  be  mounted  at  4-in.  diameter  ports 
located  at  45°  to  the  horizontal  axis  and  perpendicular  to  the  laser  optical 
axis.  Test  samples  are  introduced  through  a  4-in.  diameter  port  in  the  rear 
of  the  chamber.  This  port  contains  a  valve,  and  seal  such  that  the  samples  can 
be  introduced  with  Che  chamber  at  reduced  pressure* 

Incident  CO,,  laser  radiation,  not  falling  directly  on  the  sample,  is 
absorbed  by  a  7-in.  diameter  carbon  plate.  This  carbon  plate  is  mounted  in  a 
12-in.  diameter  aluminum  plate  which  is  supported  from  the  inner  walls  of  the 
chamber  by  a  spider. 

The  chamber  can  be  evacuated  by  diffusion  and  backing  pumps  through  a 
6-in.  diameter  port  in  the  bottom  of  the  chamber.  Dry  nitrogen  or  inert  gases 
can  be  introduced  through  an  8-in.  diameter  port  located  in  the  top  of  the 
chamber . 

2  .  CO^Laser  Transfer  Optics 

Tests  were  conducted  using  The  Aerospace  Corporation  Materials  Sciences 
Laboratory's  1-kW  CO,,  laser  which  delivered  a  multi-mode  3-in.  diameter  beam 
at  the  facility.  Alignment  of  the  optical  axis  of  the  laser  with  the  optical 
axis  of  the  test  chamber  was  accomplished  through  a  periscope.  The  periscope 
mirrors  were  8-in.  diameter  plate  copper. 


The  laser  beam  was  focused  to  a  diameter  of  V3  mm  aC  the  test  sample 
by  a  4-in.  diameter  KC&  lens  of  12-in.  focal  length.  Lens  protection  was 
accomplished  by  using  a  replaceable  3-in.  diameter  KC it  window.  Vacuum 
integrity  was  accomplished  by  O-ring  sealing  on  an  AK-coated  ZnSe  window  3  in. 
in  diameter. 


3. 


Diagnostic  s 

The  diagnostic  optical  schematic  is  shown  in  Fig.  12.  A  blackbody 
continuum  was  produced  from  a  Nernst  glow-rod  heater  (Optitron  Model  IRR-30) 
operating  at  approximately  2100°K.  The  source  was  focused  onto  a  25-mm  x 
2.5-min  slit,  through  a  Rofin  MK  II  Type  7500  light  chopper  by  a  3-in.  diameter, 
10-in.  focal  length  parabolic  mirror.  The  wavelength  region  of  4-12  pm  was 
scanned  every  1.5  sec  by  a  circular  variable  interference  filter. 

The  resulting  radiation  was  focused  as  a  12.5-mm  x  1.2-mra  slit-image  ac 
the  center  of  the  chamber  by  a  2-in.  diameter,  8-in.  focal  length  ZnSe  lens. 
Radiation  transmitted  through  the  plume  region  was  then  focused  onto  a  HgCdTe 
detector  by  a  2-in.  diameter,  2-in.  focal  length  ZnSe  lens.  The  vacuum 
integrity  as  well  as  lens  protection  were  accomplished  by  using  replaceable 
3-in.  diameter  KC&  windows. 

Data  were  processed  through  a  Princeton  HR-8  lock-in  amplifier  and 
Brush  Mark  280  strip-chart  recorder.  A  blocking  filter  with  an  8-pn  cutoff 
was  used  at  the  detector  to  prevent  spectral-order  overlap  and  to  suppress 
scattered  10.6-pm  laser  radiation  while  testing  occurred  in  the  4-8  pm 
region. 


4  .  Operating  Procedure 

Kach  experiment  was  performed  by  installing  a  carbon  rod,  evacuating 
and  refilling  (if  desired)  the  chamber,  and  then  pre-focusing  and  aligning  the 
1-kW  CO^  laser  using  a  coaligned  He-Ne  beam.  The  wavelength  range  was 
selected  and  the  appropriate  blocking  filter  was  installed.  The  circular 
variable  filter  was  rotated  at  approximately  0.7  rps.  This  typically  provided 
a  0.75-sec  4-8  pm  spectrum,  followed  by  a  filtered  7-12  pm  spectrum. 

Thus,  a  useful  spectrum  was  recorded  every  1.5  sec.  When  7-12  pm 
spectroscopy  was  selected,  the  7-12  pm  blocking  filter  acted  to  completely 
block  the  4-8  pm  data.  The  spectroscopy  was  recorded  for  a  period  of  5-10 
sec  before  the  high  power  laser  was  turned  on;  it  was  terminated  when  the 
post-burn  spectrum  resembled  the  pre-burn  spectrum  in  peak  transmission. 
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EXPERIMENTAL  RESULTS 


Our  experimental  approach  was  to  use  the  limited  time  available  with 
the  1-kW  C0?  laser  facility  to  conduct  a  survey  of  the  primary  variables 
involved . 

The  following  variables  were  investigated: 

o  Infrared  Wavelength  Range.  Two  CVIFs  were  available  for  these 
initial  experiments:  a  4-8  pm  filter  and  a  7-12  pm  filter 

o  Type  of  Carbon.  Three  carbons  were  vaporized:  Graphnol, 
pyrolytic  graphite  (perpendicular  to  lamellar  planes),  and 
Poco-graphite 

°  Probe  Spatial  Location.  1.5-sec  spectroscopy  across  the  plume  was 
performed  at  three  axial  distances  from  the  vaporizing  carbon  rod 
source:  6  mm,  18  mm,  and  30  ram 

o  Chamber  Pressure.  Experiments  were  conducted  at  pressures  of 
-2 

10  ,  1 ,  30 ,  300,  and  760  torr.  Argon  was  used  as  the  inert  gas 

A  summary  of  the  experimental  matrix  and  result  notations  appears  in 
Table  1. 

1 .  First  Experiments 

Initial  experiments  were  conducted  using  Graphnol  rods  in  a  near-vacuum 

_2 

(  10  torr  argon)  environment  and  with  the  4-8  pm  CVIF  placed  at  the 
chopper  (see  Fig.  12).  Spectral  scans  were  recorded  from  the  slit  probe-beam 
which  passed  perpendicular  to  the  incoming  10.6  Un  laser  energy,  approxi¬ 
mately  18  mm  (6  rod  diameters)  from  the  end  of  the  rod.  This  geometry  is 
shown  pictorially  shown  in  Fig.  15  and  schematically  in  Fig.  16.  As  with  most 
survey  experiments,  there  are  aLways  surprises.  Our  very  first  experiment 
(22-1)  produced  results  which  we  had  not  anticipated.  When  the  laser 

was  turned  on,  delivering  approximately  960  W  onto  the  end  of  the  3-mm 
diameter  rod,  the  4-8  pm  spectral  signature  was  quickly  almost  totally 


Table  1 


Spectroscopy  Runs 


Ran 

No. 

Chamber 

Pressure 

(torr) 

Material 

Pil cer 

x/d 

Laser 

Power 

(W) 

Comments 

Signal  Recovery 
After  Laser-Off 
( sec ) 

■ 

Oraphnol 

4-8  Um 

0  chopper 

6 

960 

10.6  um  scattering,  10-152 
sccenuacion  afcer  laser-o££ 

3 

■ 

10' 2 

Craphnol 

7-12  Um 

0  chopper 

6 

980 

10.6  um  scaccering,  10-15* 
atcenuacion  afcer  laser-off 

10+ 

2J-2 

5xlO-2 

Oraphnol 

|  7-12  um 

9  chopper 

2 

950 

Strong  scattering,  802 
attenuation  after  laser-off 

35+ 

2  +  -1 

>xl0'2 

Oraphnol 

4-8  um 

0  detector 

2 

850 

Low  scarcer,  50*  attenuacion 
during  laser-on 

30 

24-2 

10* 1 

Craphnol 

4-8  urn 

0  detector 

10 

1000 

Low  scatcer,  20*  accenuacion 
during  laser-on 

3-4 

25-1 

10'2 

Pyrolyc  ic 
Oraphice 

44-8  um 

9  detector 

2 

900 

Low  scatcer,  50*  accenuacion 
during  laser-on 

90 

25-2 

10' 1 

Pyrolyc ic 
Oraphice 

4-8  Um 

9  detector 

10 

925 

Low  scatter,  10Z  attenuation 
during  laser-on 

<1 

26-1 

5xl0”2 

Poco 

4-8  um 
@  dececcor 

2 

875 

Low  acaccer,  80-90*  attenua¬ 
tion  during  laser-on 

50  -  60 

26-2 

10'2 

Poco 

4-8  lira 

0  detector 

10 

925 

Low  scatcer,  20*  attenuation 
during  laser-on 

12 

27-1 

1 

Oraphnol 

4-8  um 
i.3  dececcor 

10 

850 

Almost  no  modulation.  10-20* 
accenuacion 

2-3 

27-2 

3<J 

Oraphnol 

4-8  Um 

0  detector 

10 

900 

Slight  modulation,  10-20* 
actenuation 

2  -  3 

27-3 

30  C 

Oraphnol 

4-8  Um 

3  dececcor 

10 

925 

Moderate  modulation,  10-20X 
at tenuation 

4-5 

28-1 

760 

Graphno l 

4-8  in 

0  detector 

6 

880 

20-40*  Absorption,  lessens 
for  t  >  3  sec 

3-4 

28-2 

760 

Pyrolytic 

Graphite 

4-8  in 

0  detector 

6 

980 

Modulation  (scattering), 
spectrum  recovery  at  ~3  aec 

3 

29-1 

760 

Poco 

4-8  um 
(p  dececcor 

6 

Severe  modulation  (scatter¬ 
ing),  spectrum  reappears 
t  >  3  sec 

6 
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Fig.  16.  Probe  Bean  Geometry 


masked  by  particle-modulated  10.6  im  scattered  light  from  the  CC>2  laser. 
Unexpectedly,  there  were  sufficient  carbon  particles  near  Che  end  of  the 
carbon  rod  source  to  produce  the  amount  of  spectral  modulation  that  was 
observed.  What  had  been  anticipated  was  that  carbon  particle  conden  ation 
would  occur  some  distance  from  the  carbon  source.  Based  upon  this  assumption, 
we  expected  that  there  would  be  a  condensation  zone  which  would  cover  the 
range  from  gaseous  carbon  through  small  particles  to  a  coalesced  particle 
regime. 

With  our  geometrically  well-defined  optical  diagnostic  capability,  if  a 
vapor  zone  had  existed  in  this  region,  there  would  have  been  little  attenuation 
of  che  transmitted  spectrum  and  a  low  level  of  observed  10.6  pm  scattered 
light.  Our  observations  are  presented  in  Fig.  17  which  shows  a  typical 
pre-burn  spectrum,  a  portion  of  che  burn  spectrum,  and  part  the  first 
post-burn  spectrum  which  demonstrates  the  attenuation  of  the  transmitted  beam 
tiirough  the  cooling  plume. 

From  this  first  test  result,  we  concluded  that  there  were  a  substantial 
number  of  particulates  18  mm  from  the  end  of  the  carbon  source.  Also,  the 
intensity  of  the  scattered  light  was  strong  from  the  point  of  turn-on  and 
remained  fairly  conscant  through  the  4-sec  burn  sequence.  In  the  first 
post-burn  spectrum,  where  the  10.6-pra  scattered  radiation  was  no  longer 
present,  a  10  to  15  percent  attenuation  was  observed  at  all  wavelengths. 
Complete  return  to  the  pre-burn  spectrum  occurred  in  about  3  sec.  The  peak 
transmission  points  on  each  side  of  the  atmospheric  6.25-pm  water  absorption 
correspond  to  4.86  pro  and  6.86  pm.  Approximately  10  percent  greater 
absorption  was  observed  from  the  longer  wavelength  in  the  first  post-burn 
s  pec  c  rum. 

A  second  experiment  (23-1)  was  conducted  under  the  same  conditions, 
using  the  7-17.  pm  CV1F.  All  the  observations  noted  in  the  first  experiment 
apply  also  to  the  results  of  this  experiment,  including  the  increased 
attenuation  at  longer  wavelengths.  In  both  experiments,  there  were  no  intense 
or  easily  identified  molecular  or  structured  absorption  features  in  the 
post-burn  spectra. 


Spectroscopy 


The  last  experiment  in  this  initial  series  (23-2)  was  a  repeat  of  23-1 
with  diagnostics  conducted  6  mm  (compared  with  18  mm)  from  the  end  of  the 
carbon  source.  Representative  7-12  tin  spectral  scans  from  this  run  are 
shown  in  Fig.  18.  It  was  expected  that  there  would  be  fewer  particles  at  this 
distance  and  therefore  less  absorption  and  less  scattering  of  the  10.6  tin 
C02  laser  energy. 

The  exact  opposite  was  observed.  There  was  greater  attenuation  of  the 
diagnostic  beam — approximately  80  percent  for  the  first  post-burn  spectral 
scan.  Tne  intensity  of  the  modulated  10.6-tnn  light  was  50  percent  greater. 
Again,  there  was  stronger  attenuation  at  long  wavelengths.  Measured 
attenuation  ratios  were  1.72  between  7  and  12  pn.  This  number  compares  well 
with  a  ratio  of  1.71  if  the  attenuation  is  proportional  to  wavelength. 
Surprisingly,  the  recovery  to  the  pre-burn  transmission  conditions  required 
over  35  sec  indicating  that  particulates  were  either  continuing  to  be 
generated  after  the  00^  laser  was  shut  off  or  required  long  settling  times 
to  clear  the  field  of  view. 

On  the  assumption  that  vacuum  plume  streamlines  near  the  end  of  the 
vaporizing  carbon  rod  "loft"  small  particles  which  settle  when  the  laser  is 
turned  off,  settling  times  were  calculated  based  upon  both  continuum  and 
collisionless  gas  dynamics.  Details  of  the  calculations  are  given  in  Appendix 
B.  The  conclusion  is  that  if  the  long  post-burn  spectral  recovery  times  are 
caused  by  the  slow  settling  of  lofted  particles,  Chen  the  particles  arc 
exceedingly  small  (10-30  nm) . 

After  che  results  of  the  first  three  tests  were  assessed,  a  new 

experimental  strategy  was  developed.  It  was  decided  that  our  best  chance  of 

reducing  the  scattered  10.6  pm  lighc  and  performing  spectroscopy  during  the 

burn  phase  would  be  to  perform  4-8  spectroscopy  on  the  remaining 

experiments.  The  4-8  |jn  blocking  filter,  which  highly  attenuates  10.6  pm 

light,  was  moved  to  a  position  directly  in  front  of  the  detector.  Using  this 

strategy,  we  completed  a  series  of  12  additional  experiments  to  study 

different  types  of  graphite  with  near  and  far  field  diagnostics  at  ambient 

-2 

pressures  ranging  from  10  to  760  torr  (see  Table  l).  These  experiments 
are  described  in  the  following  subsections. 
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4. 


Spectroscopy 


Effect  of  Varying  Graphite  Material  and  Probe  Axial  Location 

In  this  series,  six  experiments  (24-1  through  26-2)  were  conducted  at 
near-vacuum  chamber  conditions.  Three  graphites— Graphnol  N3M^,  pyrolytic 
graphite  (Super  Temp,  continuously  nucleated),  and  Poco  (AXF.Q1) — were  studied 
at  diagnostic  probe  points  6  and  30  mm  from  the  end  of  the  carbon  rod.  The 
results  are  partly  predictable  and  otherwise  interesting.  The  10.6  urn 
scattered  radiation  problem  was  substantially  reduced,  allowing  4~8  pm 
spectral  scans  to  be  seen  during  the  laser-on  time.  The  data  for  the  two 
diagnostic  probe  distances  confirmed  our  earlier  conclusion;  namely,  that 
greater  numbers  of  particles  are  found  closer  to  the  end  of  the  plume  source. 
Also,  as  in  the  initial  experiments,  the  post-burn  spectral  recovery  period  is 
much  longer  during  probing  close  to  the  rod,  i.e. ,  3-12  sec  at  a  distance  of 
30  mm  and  30-90  sec  at  6  mm.  The  percentage  of  probe  energy  attenuacion  is  3 
to  20  percent  at  30  mm;  at  6  mm  from  the  end  of  the  rod,  as  much  as  80  to  90 
percent  attenuation  was  observed. 

In  many  ways,  the  three  types  of  graphite  tested  gave  similar  results. 
Absorption  and  scattering  of  infrared  radiation  by  particles  was  observed  for 
all  three  materials.  At  an  axial  station  of  30  mm  (10  rod  diameters),  the 
fractional  absorption  (up  to  20  percent)  was  essentially  independent  of 
material  type  (see  Fig.  19).  This  figure  also  illustrates  that  during  the 
run,  the  level  of  modulated  10.6-pm  light  is  least  for  Graphnol  and  greatest 
for  pyrolytic  graphite. 

At  6  mm  (see  Fig.  20),  much  greater  absorption  of  the  probe  energy 
occurred.  All  three  graphites  showed  an  increased  particle  density  build-up 
with  time.  The  build-up  was  fastest  for  pyrolytic  graphite  because  of  its  low 
thermal  conductivity  along  the  rod  axis  and  the  resulting  much  shorter  time 
(V).l  sec)  to  steady-state  ablation.  The  greatest  absorption  during 
laser-on  (^90  percent)  occurred  with  the  Poco-graphite . 
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Fig.  20.  Laser-on  Spectroscopy  Close  to  End  of  Rod 


Post-burn  spectral  signal  recovery  times  were  approximately  30,  60,  and 
90  sec  for  Graphnol,  Poco,  and  pyrolytic  graphite,  respectively.  The  slow 
signal  recoveries  can  be  seen  in  the  first  three  post-burn  spectral  scans  for 
these  materials  shown  in  Fig.  21.  We  note,  cautiously,  chat  all  three 
graphites  indicate  the  presence  of  an  absorption  band  at  h.7.4  pra  (1350 
cm  *■)  in  the  immediate  post-burn  spectra.  There  has  been  no  attempt  to 
assign  this  feature  to  a  molecular  species. 

Overall,  these  experiments  indicate  that  more  particulates  (absorption) 
were  generated  with  the  Poco  graphite  than  with  the  other  two  graphites. 
However,  the  modulated  10.6-pm  light  (scattering)  was  greatest  for  the 
pyrolytic  graphite. 


3.  Effect  of  Finite  Ambient  Pressures 


A  series  of  three  experiments  (27-1,2,3)  was  run  with  Graphnol  rods  at 
argon  pressures  of  1,  30,  and  300  torr.  A  diagnostic  distance  of  30  mm  was 
chosen,  since  an  increase  in  particulate  concentration  was  anticipated  due  to 
the  higher  pressures  and  mixing  of  the  hot  carbon  vapors  with  the  cold  ambient 
gas. 


Once  again,  the  results  of  these  tests  did  not  confirm  our  expectations. 
Atcenuation  of  the  probe  beam  was  approximately  20  percent  during  the  laser-on 
time  and  was  essentially  independent  of  pressure.  However,  modulation  of  the 
10.6-pm  light  was  ^3.5  times  larger  at  300  torr  chan  at  1  torr.  These 
effects  are  illustrated  in  Fig.  22.  There  were  indications  that  the  7.4-pm 
band  was  again  present  and  particularly  strong  in  some  of  Che  laser-on 
spectra.  There  may  also  be  a  spectral  feature  at  5  pm  which  correlates  with 
the  7.4-pm  band  intensity.  The  triatomic  carbon  molecule,  C  ,  has  a 
rocat ion/vibration  band  near  5  pm.  Spectral  recovery  times  were  on  the 
order  of  3  to  5  sec. 

A  final  series  of  three  experiments  (28-1,  28-2,  29-1),  using  all  three 
graphites,  was  carried  out  with  an  ambient  pressure  of  760  torr  in  argon.  The 
4-8  pm  spectroscopy  was  erformed  6  ram  from  the  end  of  the  carbon  source. 
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SPtCTRUM  POST-BURN 


Transmission  Recovery  Effects 


This  series  showed  interesting  results.  First,  the  Graphnol  produced  the 
least  amount  of  modulation  of  10.6  urn  radiatio:..  The  attenuation  increases 
with  time  for  about  3  sec  and  then  decreases.  Poco-graphite  produced  the 
greatest  10.6  Mm  modulation.  However,  after  about  3  sec  of  burn,  the 
modulation  decreased  and  the  spectrum  began  to  reappear.  Maximum  attenuation 
in  the  post-burn  spectra  appeared  to  be  20  to  30  percent.  Pyrolytic  graphite 
showed  a  similar  effect  with  a  spectrum  reappearing  approximately  2  sec  into 
the  burn.  These  effects  are  shown  in  Fig.  23. 

These  tests  apparently  indicate  that  about  halfway  into  the  burn,  all 
the  graphites  produced  sufficient  particulate  densities  to  initiate  blockage 
of  the  incoming  CO^  laser  energy.  Many  large-scale  (0.1-1  nun)  "fluffy" 
particles  of  carbon  condensate  were  observed  in  the  chamber  after  each  of 
these  runs.  Some  of  these  were  collected  for  later  analysis.  The  post-test 
appearance  of  the  samples  (Fig.  24)  indicates  that  the  major  source  of  these 
large  pieces  may  have  been  condensation  and  subsequent  detachment  from  the  end 
of  the  vaporizing  (subliming)  carbon  rod  during  the  run.  Graphite  rods 
exposed  to  laser  radiation  at  lower  ambient  chamber  pressures  did  nbt  exhibit 
this  phenomenon. 

D.  OBSERVATIONS 

At  least  three  unexpected  results  were  produced  from  this  series  of 
experiments : 

o  Attenuation  of  the  probe  beam  by  the  plume  increased  as  the  beam  was 
positioned  closer  to  the  end  of  the  carbon  rod. 

o  Transmission  measurements  through  the  plume  showed  preferentially 
increased  attenuation  at  longer  wavelengths. 

o  Extremely  long  (up  to  90  sec)  probe  beam  attenuation  recovery  times 
were  observed  after  laser  turn-off. 

Understanding  these  results  is  an  important  goal,  since  they  appear  to  be 
related  to  each  other  and  should  reveal  a  great  deal  about  the  fundamental 
processes  involved. 


Spectroscf-  ,  Chamber  Pressure 


Stronger  attenuation  of  the  probe  beam,  when  it  was  positioned  closer 

to  the  end  of  the  carbon  rod,  can  only  be  interpreted  as  more  particles  per 
2 

cm  .  Thus,  the  data  indicate  that  if  a  condensation  zone  exists,  it  is 
within  the  6mm  minimum  distance  studied  in  our  experiments.  The  reason  for 
this  is  depicted  qualitatively  in  Fig.  25.  The  upper  figure,  25a,  shows  the 
behavior  of  the  driving  force  for  condensation  processes  (i.e.,  the 
supersaturation  ratio)  as  a  function  of  distance  from  the  end  of  the  rod. 
Figure  25b  illustrates  Che  behavior  of  the  collision  frequency  of  molecules 
within  the  plume.  Since  the  number  of  condensed  particles  present  at  any  time 
depends  on  the  product  of  these  two  quantities,  the  particle  density  reaches  a 
maximum  value  and  then  decreases  due  to  volumetric  dilution  (Fig.  25c).  If 
condensation  is  occurring  in  our  experiments,  the  peak  particle  density  is 
very  close  to  the  end  of  the  rod. 

A  very  reasonable  and  likely  alternative  explanation  of  what  we  are 

18  19 

seeing  is  that  the  particles  are  being  emitted  ’  from  the  vaporizing 
carbon  surface.  In  this  instance,  we  would  expect  to  see  the  attenuation 
continue  to  increase  as  the  probe  beam  approaches  the  end  of  Che  rod.  This  is 
shown  schematically  by  the  dashed  line  in  Fig.  25c.  More  detailed  probing 
with  greater  resolution  near  the  end  of  Che  vaporizing  rod,  some  carefully 
planned  experiments,  and  particle  collection  will  be  required  to  finally 
resolve  the  question  of  particle  origins. 

The  observed  increased  attenuation  with  wavelength  and  the  long  signal 
recovery  times  after  laser-off  are  both  connected  in  some  way  with  particle 
sizes.  Our  calculations  of  particle  sizes,  based  upon  the  assumption  that 
they  are  "lofted"  within  the  chamber  (Appendix  B) ,  indicate  particle  sizes  the 
order  of  0.001  times  the  wavelength  of  the  probe  light.  In  this  regime,  we 


1  Whittaker ,  A.  G.  and  P.  L.  Kintner,  "Particle  Emission  and  Related  Morpho¬ 
logical  Changes  Occurring  During  the  Sublimation  of  Graphitic  Carbons," 
Carbon,  Vol.  14,  1976,  pp.  257-265. 

l^Abrahatnson f  j.f  "Graphite  Sublimation  Temperatures,  Carbon  Arcs,  and 
Crystallite  Erosion,"  Carbon,  Vol.  12,  1974,  pp.  111-141. 
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would  expect  to  see  absorption  scale  inversely  with  wavelength-  rather  than 

directly  with  wavelength  as  observed  experimentally.  However,  if  the  particles 

were  slightly  larger  than  the  wavelength  scanned,  the  observed  wavelength 

2 1 

dependence  might  be  explainable  in  terms  of  the  Mie  theory.  Production  of 
large  particles  is  consistent  with  the  physical  properties  of  the  graphites 
used.  Representative  grain  sizes  for  Poco  graphite,  Graphnol,  and  pyrolytic 
graphite  are  the  order  of  3ym,  30tfn,  and  300ym,  respectively.* 


Finally,  we  note  an  unresolved  discrepancy  in  the  above  discussion. 

If  the  particles  are  large  so  as  to  possibly  explain  the  attenuation  behavior 
with  wavelength,  they  will  quickly  fall  to  the  bottom  of  the  chamber  after 
laser-off.  However,  very  small  particles  are  required  to  explain  the  long 
signal  recovery  times.  Alternative  explanations  are  needed. 


We  plan  to  resolve  these  qustions  in  the  near  future  by  carrying  out 
additional  experiments  using  the  hardware  and  diagnostic  methods  now 
available.  The  primary  objectives  of  new  work  would  be  to: 


o  Identify  the  particle-generating  mechanism  (emitted,  condensed,  or 
both) 

o  Collect  and  characterize  particles  over  a  wide  range  of  ablating 
surface  temperatures 

o  Relate  observed  particle  effects  (i.e.,  absorption/scattering)  to 
the  solid-state  properties  of  the  type  of  carbon  vaporized  and  to 
the  properties  of  collected  particles 

o  Further  investigate  possible  molecular  absorption  features  seen 
near  5ym  (likely  triatomic  carbon,  C3)  and  7.4  ym 

o  Extend  our  diagnostics  to  cover  the  entire  2.5-14  pra  raid-IR 
spectra  1  range 


^Hottel,  h.  c.  and  A.  F.  Sarofin,  Radiative  Transfer,  McGraw-Hill  Book  Co., 
N.Y.,  1967,  Chapter  12,  13. 

2  Wan  der  Hulst,  H.  C. ,  Light  Scattering  by  Small  Particles,  Dover  Publi¬ 
cations,  Inc.,  NY  ( 1981) . 

*  Chase,  A.B.,  Personal  communication,  The  Aerospace  Corporation,  El  Segundo, 
CA. 
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V.  SUMMARY  AND  CONCLUSIONS 


The  overall  objectives  of  oar  four  years  of  work  have  been  to  gain  a 
better  understanding  of  the  superior  energy-dissipating  efficiency  of  carbon 
materials  and  to  determine  the  limit3  of  their  outstanding  thermal  protection 
system  performance.  The  method  used  to  accomplish  this  goal  has  been  to 
develop  analytical  techniques  which  allow  detailed  analyses  of  nonequilibrium 
carbon  laser  ablation  experiments.  Our  analytical  approach,  data  analysis 
conclusions,  and  their  relevance,  are  summarized  below. 

In  order  to  relate  experimentally  measured  quantities,  such  as  mass 
loss  rates  and  surface  temperatures,  to  desired  thermochemical  properties,  we 
developed  a  multispecies  nonlinear  model  of  nonequil ibriura  carbon  ablation 
under  laser  irradiation.  The  initial  model  was  based  upon  solving  the  species 
mass,  momentum,  and  energy  conservation  equations  across  the  nonequilibrium 
region  (Knudsen  layer)  at  the  surface.  The  importance  of  ad  hoc  closure 
methods  associated  with  the  conservation  equation  approach  was  investigated  by 
using  higher  order  methods  (Boltzmann  equation  and  Molecular  Dynamics)  to 
solve  for  the  changes  in  macroscopic  variables  across  the  Knudsen  layer. 

These  methods  established  the  validity  of  the  more  approximate  conservation 
equation  solutions  which  were  used  for  our  data  analyses.  We  have  also  now 
established  a  hierarchy  of  models  which  can  be  applied  to  other  complex  and 
theoretically  difficult  physical/cheraical  problems. 

Using  the  above  model  as  a  data  analysis  tool,  we  find  that  our 
continuous  wave  (CW)  and  pulsed  laser  ablation  experiments  have  indicated  the 
following: 

o  The  vapor  pressure  of  carbon  is  best  represented  by  the  Lawrence 
Livermore  Laboratory  thermochemical  properties . ^ ^ 

o  The  melt  temperature  of  carbon  is  in  the  range  of  4500-4800°K 

o  Vaporization  coefficients,  which  define  the  kinetic  races  of 
carbon  vaporization,  are  strongly  temperature-dependent.  The 
vaporization  coefficients  for  Graphnol  and  for  pyrolytic  graphite 
(when  vaporization  takes  place  perpendicular  to  the  planes)  are 
about  the  same.  When  pyrolytic  graphite  vaporization  occurs 
parallel  to  the  planes,  the  vaporization  coefficients  are  approxi¬ 
mately  a  factor  of  three  higher. 


o 


The  high  temperature  ablated-surface  eruissivities  of  Graphnol  and 
pyrolytic  graphite  (viewed  parallel  to  the  planes)  are  near 
unity.  For  pyrolytic  graphite  viewed  perpendicular  to  the  planes, 
the  high  temperature  ablated-surface  eraissivity  is  estimated  from 
our  data  to  be  0.85. 

Our  analytical  model,  together  with  the  above  carbon  properties  data, 
allows  reliable  thermal  protection  system  design  calculations  to  be  made  for 
the  projected  use  of  carbon  materials  in  very  severe  heating  and 
none qui librium  ablation  environments.  The  most  serious  limitations  to 
ultimate  performance  capability  appear  to  be:  degradation  of  mechanical 
properties  at  temperatures  above  400Q°K,  rapidly  increasing  surface  recession 
rates  which  increase  from  25  mm/sec  (1  in/sec)  to  300  mm/ sec  (1  ft/sec)  as  the 
cemperature  rises  from  4200  to  4500°K,  and  melting  above  4500°K.  To  reach 
450O°K,  a  local  heat  flux  of  about  2  MW/cm^  (1,760,000  Btu/ft^-sec)  is 
required. 

In  addition  to  the  above  basic  studies,  a  preliminary  investigation  of 
particle  effects  in  a  laser-vaporized  carbon  plume  has  been  completed  using 
spectroscopic  methods.  From  these  initial  experiments,  the  following 
observations  have  been  made: 


o  Severe  particulates  blockage  (80  to  90  percent)  of  probe-radiation 
from  an  infrared  source  occurs  when  one  views  through  the 
laser-vaporized  carbon  vapor  plume  at  a  location  close  (6  mm)  to 
the  vaporizing  surface  (p  **  10-2  torr). 

CO 

o  The  absorption  increases  with  wavelength  and,  in  some  instances, 
is  proportional  to  wavelength. 

c  When  probes  are  made  at  locations  farther  from  the  source  (18  ram, 
30  mm),  the  blockage  is  much  less. 

o  During  the  laser  heating  time,  the  spectroscopic  signature  is 
strongly  modulated  by  10.6  laser  radiation  scattered  by 
particles  in  the  plume. 

The  likely  source  of  the  particles  is  particulate  emission  from 
the  ablating  surface.  Condensation  in  the  highly  supersaturated 
carbon  vapor  remains  as  a  possible  alternate  or  simultaneous 
particle  generating  source. 


o 


Although  quantitative  comparisons  are  difficult,  Poco  graphite 
appears  to  produce  more  particles  (and  greater  infrared  radiation 
blockage)  than  Graphnol  or  pyrolytic  graphite. 

Our  present  experiments  and  analysis  give  inconsistent  conclusions 
regarding  particle  sizes. 

At  higher  chamber  pressures,  especially  at  1  atra,  more  particles 
and  greater  laser  radiation  blockage  were  observed. 
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APPENDIX  A 


ONE-DIMENSIONAL  STEADI-STATE  ABLATION 
OF  A  SLENDER  ROD 


The  components  of  the  steady-state  energy  balance  for  laser  irradiation 
of  an  absorbing  solid  material  target  are  shown  in  Fig.  A-l.  In  order  to  gain 
a  complete  understanding  of  the  energy  dissipation  efficiency  of  a  particular 
material,  one  must  predict  and  experimentally  measure  each  of  these 
components.  In  many  instances,  che  conduction  of  energy  into  the  solid, 

^cond’  c*oe9  noc  achieve  a  steady-state  (constant)  value  during  the  time  tnat 
the  materiaL  is  irradiated  by  the  laser.  This  complicates  the  analysis 
process  by  requiring  that  predicted  quantities  and  experimental  data  must  both 
Do  determined  as  functions  of  time. 

Often  the  transient  heating  time  is  lengthened  by  radial  conduction  of 
energy  away  from  a  small  laser- irrad iated  spot  on  the  surface.  In  9uch 
instances,  at.  .  ate  analysis  is  further  complicated  by  the  laser-ablated  hole 
which  forms  in  the  solid  and  which  may  alter  the  expansion  of  vaporized 
material  away  from  the  surface.  A  slender  rod  sample-geometry  for  laser 
irradiation  experiments  eliminates  radial  conduction  and  hole-boring  effects. 
The  end  of  the  rod  is  "flood-loaded"  oy  the  incoming  laser  energy,  and  a 
steady-state  ablation  rate  and  temperature  distribution  are  rapidly 
uutabl  ished. 

T;ie  purpose  of  this  appendix  is  to  describe  the  differential  equation 
and  boundary  conditions  governing  this  problem,  to  show  numerically  obtained 
representative  solution  behavior,  and  to  briefly  discuss  two  applications  of 
thi^  type  of  solution. 

DIFFERENTIAL  EQUATION  AND  BOUNDARY  CONDITIONS 

The  he.it  conduction  equation  describing  energy  transport  for  a 
one-dimensional  steady-state  ablation  process  was  first  derived  by 


Fig.  A-l.  Energy  Partitioning  at  Surface 
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Landau. ^  For  stationary  surface  recession  velocity  s,  and  distance  x 
measured  from  the  moving  surface,  Che  governing  differential  equation  is 


(  A —  1 ) 


where  (Xp  is  the  thermal  diffusivity  and  T  is  Che  temperature.  Landau 
presented  solutions  to  this  equation  with  appropriate  boundary  conditions.  In 
Ref.  12,  Landau's  solution  was  extended  to  include  a  melted  layer  of  absorbing 
and  conducting  liquid  above  the  solid. 

For  a  slender  cylindrical  rod,  the  one-dimensional  heat  conduction 
equation  when  similarly  transformed  becomes 


4ocr 


(  A-2  ) 


wnere 


=  thermal  conductivity 
=  rod  diameter 
=  emissivity  of  rod  surface 

■  Stefan-Boltzmann  constant  for  radiative  heat  transfer 


The  term  on  the  right-hand  side  of  Eq.  (A-3)  arises  due  to  radiative  loss  of 
energy  from  the  cylindrical  surface  of  the  rod. 


When  we  non-dimensionalize  by  the  end-surface  temperature  Tg  and  the 
thickness  of  the  thermal  layer  (5,  it  is  convenient  to  define  new  variables  T 
and  x  by 


T  =  I- 


(A-J) 


Equation  (A-3)  now  becomes 
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subject  to  the  boundary  conditions 

T(0)  =  1  ,  T( 1)  =  T^ 


(A-4) 


(A-5) 


Because  the  thickness  of  the  thermal  layer  6  appears  directly  in  the 
differential  equation,  a  further  condition  is  needed  to  define  6.  A  global 
energy  balance  provides  the  required  constraining  equation,  i.e. 
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The  term  on  the  Left  in  Eq.  (A-6)  is  Che  energy  conducted  into  the  rod  at  the 
end.  This  energy  is  dissipated  in  two  ways.  The  integral  term  on  the  right 
represents  the  total  energy  radiated  from  the  cylindrical  surface  of  the  rod. 
The  last  term  on  the  right-hand  side  is  the  energy  required  to  heat  the  rod 
from  temperature  T^  to  Tg ,  assuming  that  Cp  is  not  a  function  of 
temperature . 


A -2  .  NUMERICAL  SOLUTION  PROCEDURE 

Equation  (A-4)  is  a  nonlinear  two-point  boundary  value  problem  with 

boundary  conditions  [Eq.  (A-5)].  In  the  present  analysis  it  was  solved  using 

standard  second  order  accurate  centered  differencing  after  the  application  of 

quas  i  1  inear izat  ion  to  obtain  a  linear  approximation.  For  any  fixed  5  >  0, 

the  iterations  associated  with  the  quasilinearization  converge  fairly 

—4 

rapLdly.  Iteration  tolerances  of  10  were  employed  with  a  max-norm  con¬ 
vergence  criterion.  Typically,  between  two  and  six  iterations  were  needed. 


The  grid  spacing,  h,  used  in  most  runs  was  h  =  0.02.  A  few  runs  were  made 
with  h  =  0.01  and  h  =  0.005  to  verify  grid  function  convergence.  IL  was  found 
that  the  h  =  0.02  results  were  within  a  few  percent  of  converged  solutions. 

The  energy  balance,  Eq.  (A-6),  provides  the  equation  for  determining  6. 
This  expression  is  highly  nonlinear  with  respect  to  6,  because  T  and  dT/dx 
are  implicitly  functions  of  6  through  Eq .  (A-4),  in  addition  to  the  explicit 
quadratic  dependence  shown  in  Eq.  (A-6).  Equation  (A-6)  was  solved  by 
Newton's  method,  but  with  numerical  approximations  of  derivatives  with  respect 
to  6  used  for  terms  for  which  analytical  derivatives  could  not  be  calculated 
(the  terms  containing  T  and  dT/dx).  Because  of  Che  need  for  such  approxi¬ 
mations,  and  also  because  of  the  high  degree  of  nonlinearity,  convergence  was 
sometimes  slow.  Nevertheless,  complete  solutions  to  the  system  consisting  of 
Eqs.  (A-4),  ( A— 5 ) ,  and  (A-6)  could  always  be  obtained  within  approximately  20 
iterations  of  Eq.  (A-6). 

The  overall  computational  algorithm  consists  of  the  following  steps: 

a.  Guess  a  value  of  <S  and  a  temperature  distribution,  T(x). 

b.  Solve  Eq.  (A-4)  with  boundary  conditions  (A-5)  (Note  that  Che 
quasilinear  iterations  should  be  carried  to  convergence.). 

c.  Perform  a  single  Newton  iteration  step  on  Eq.  (A-6)  to  obtain  an 
improved  value  of  6. 

d.  Test  convergence  of  6.  If  6  is  converged,  stop.  Otherwise, 
repeat  steps  2  and  3. 


A-3. 


REPRESENTATIVE  SOLUTIONS 


Anticipating  Che  use  of  calculated  solutions  in  conjunction  with  the 
experiments  described  in  Section  V  of  this  report,  we  selected  a  nominal  rod 
diameter  D  of  3mm  (9.84  x  10  ^  ft).  The  values  of  heat  capacity  and 
thermal  conductivity  It  used  for  initial  calculations  were 

C  =  0.543  Btu/ lb-° R 

P 

k  =  0.00295  Btu/ft-sec-°R 

These  values  are  representative  of  graphite  at  high  temperatures.  In  order  to 
minimize  the  time  required  to  calculate  solutions,  the  surface  recession  rate 
s  (ft/sec)  was  represented  as  a  function  of  surface  temperature  Tg  by  the 
f uncc ion 


-E/KT 

s  =  Ae  3 

where  A  =  7.69  x  10^  and  E/K  a  187,207  °R  for  T  in  °R.  This  expression 
was  obtained  from  carbon  nonequilibrium  ablation  calculations  discussed  in 
Kefs.  13-15. 

From  the  calculated  steady-state  temperature  distributions  in  the  rod, 
Che  heat  conducted  into  the  rod  q%  and  the  thermal  layer  thickness  6  were 
determined.  Representative  values  of  q^  and  6  for  surface  emissivit ies  of 
0. 5-1.0  at  increments  of  0.1  and  for  surface  temperatures  of  5400°R  (3000°K)  , 
6000° K  (3333°K) ,  6600°K  (3667°)  and  7200°P  (4000°K)  are  shown  in  Table  A-l. 

For  the  lowest  surface  temperature,  the  surface  recession  rate  a  is 
essentially  zero.  In  this  case,  the  solution  represents  a  steady  stace  energy 
balance  in  which  laser  energy  is  conducted  along  the  rod  and  radiated  out  the 
si-Je.  For  the  highest  surface  temperature,  the  surface  recession  rate  s  is 
11.97  mm/sec  (4.71  in. /sec).  In  this  instance,  considerable  laser  energy  is 
accommodated  by  the  ablation  process,  and  the  remaining  energy  is  conducted 
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Table  A-l.  Dependence  of  Calculated  Results  on  Surface 
Emissivity  and  Temperature 


d*d0 

=  0.00984  ft 

Btu/ ft  ^se  c 

T  (°R) 

e 

&  (ft) 

• 

5400 

1.0 

0. 185781 

1012.876 

0.9 

0.194788 

960.385 

0.8 

0.205389 

904.922 

0.7 

0.218107 

845.907 

0.6 

0.233724 

782.528 

0.5 

0.253395 

713.587 

6000 

1.0 

0. 168189 

1343.513 

0.9 

0.177473 

1275.782 

0.8 

0.  188365 

1204. 191 

0.7 

0.201642 

1127.959 

0.6 

0.218171 

1046.062 

0.5 

0.239578 

957.011 

6600 

1.0 

0.250717 

2276.801 

0.9 

0.268295 

2203.255 

0.8 

0.289088 

2126.875 

0.7 

0.314296 

2047.339 

0.6 

0.345518 

1964.284 

0.5 

0.385383 

1877.251 

7200 

1.0 

0.111267 

16137.18 

0.9 

0. 112002 

16121.95 

0.8 

0. 112747 

16106.71 

0.7 

0.113502 

16091.46 

0.6 

0.114268 

16076.21 

0.5 

0.115044 

16060.95 
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into  the  rod  and  again  radiated  out  the  aide.  We  note  that  when  8  is  high, 

<S  and  are  insensitive  to  the  emissivity  e;  whereas  when  s  approaches 
zero,  6  increases  and  qc  decreases  as  £  decreases,  for  a  given  surface 
temperature  Tg.  The  sensitivity  of  calculated  solutions  to  the  rod  diameter 
and  the  therraaL  conductivity  is  shown  in  Tables  A-2  and  A-3,  respectively. 
Results  are  again  shown  for  the  same  surface  temperatures  as  in  Table  A-l.  To 
a  rough  approximation,  except  for  T^  3  7200°R,  the  dependence  on  rod 
diameter  and  thermal  conductivity  is  represented  by 


(effecc  of  rod  diameter) 


-1/2 


(effect  of  conductivity) 


1/2 


where  d  and  k  are  the  nominal  values  used  for  the  calculated  results 
o  o 

given  in  Table  A-l.  For  T  **  7200°R,  6  and  q  are  relatively 
insensitive  to  diameter;  q^  increases  slightly  as  k  increases,  however,  6 
is  nearly  directly  proportional  to  k,  varying  by  a  factor  of  70  when  k 
increases  two  orders  of  magnitude. 


A -4.  APPLICATIONS 


The  application  of  these  solutions  to  data  analysis  is  straightforward 
and  direct.  Unfortunately,  we  were  not  able  to  obtain  measured  surface 
temperatures  on  the  ablating  end  of  the  rod  during  the  test  series  described 
in  Section  V  of  this  report.  For  low  laser  energy  input  levels,  we  believe 
that  these  solutions  can  be  used  to  determine  the  emissivity  £  of  the 
graphite  rod.  For  the  higher  laser  input  energies,  this  solution  method 
provides  a  tool  which  can  be  readily  used  to  evaluate  each  term  (see  Fig.  A-l) 
in  the  energy  balance  equation  for  steady-state  ablation.  From  this,  the 
contribution  of  each  dissipation  mechanism  to  the  overall  ablation  efficiency 
can  be  determined  for  various  materials  of  interest. 


Table  A-2.  Sensitivity  to  Rod  Diameter 


£  =  0.9 

T 

d/d 

6 

• 

q 

o 

c 

5400 

0.  1 

0.064249 

3054.506 

l 

0.194788 

960.385 

2 

0.284471 

681.461 

3 

0.347652 

556.364 

10 

0.634412 

304.924 

6000 

0. 1 

0.056490 

3997.324 

1 

0. 177473 

1275.782 

2 

0.253511 

909.671 

3 

0.312051 

746.761 

10 

0.601726 

419.500 

6600 

0.1 

0.057084 

3515.368 

l 

0.268295 

2203.255 

2 

0.410115 

1831.997 

3 

0.514038 

1688.795 

10 

0.872119 

1465.880 

7200 

0.  1 

0.074029 

17339.88 

1 

0.112002 

16121.95 

2 

0. 1 L  5529 

16052.46 

3 

0.116719 

16029.60 

10 

0.118440 

15997.44 

Table  A-3.  Sensitivity  to  Thermal  Conductivity 


k/k 

0 

6 

• 

q 

c 

■ 

k/k 

o 

<5 

• 

qc 

0.  1 

0.063939 

305.535 

6600 

0.  1 

0.087279 

1465.801 

i 

0.194788 

960.385 

i 

0.268295 

2203.255 

10 

0.633841 

3043.923 

10 

0.561721 

5496.961 

0. 1 

0.062899 

423.384 

7200 

0.1 

0.011900 

15992.44 

1 

0. 177473 

1275.782 

l 

0.112002 

16121.95 

10 

0.555039 

3979.489 

—  .  .  -i 

_ I 

10 

i 

0.741124 

17338.61 

APPENDIX  B 


PARTICLE  SIZES  DERIVED  FROM  SETTLING  TIMES 

If  it  is  assumed  that  the  long  probe-beam  recovery  times  seen  experi¬ 
mentally  are  associated  with  particles  "lofted1'  in  the  chamber  at  the  probe 
station,  then  approximate  particle  sizes  can  be  obtained  by  comparing 
predicted  settling  times  with  those  seen  in  the  data.  Alternatively,  the 

settling  time  t,  the  settling  velocity  V  ,  and  the  settling  length  L  are 

s  s 

related  by  V  =  L  /t.  Whan  L  =  1  ft  (half  the  chamber  diameter)  is 
s  s  s 

used  for  the  settling  length,  the  observed  long  settling  times  of  30  to  90 
sec,  give  values  from  0.011  to  0.033  ft/sec.  The  purpose  of  this 

appendix  is  to  determine  for  what  range  of  particles  sizes  the  predicted 
velocities  are  in  the  observed  range.  Calculations  based  upon  the  bounding 
cases  of  concinuum  and  collisionless  gas  dynamics  were  carried  out. 

3-1.  PARTICLE  SIZES  PREDICTED  FROM  CONTINUUM  GAS  DYNAMICS 


For  a  particle  settling  in  a  fluid,  the  terminal  velocity  is  determined 

by  a  balance  of  the  acceleration  force  due  to  gravity  and  the  retarding 

force  F^  due  to  viscous  forces.  For  a  particle  of  radius  R  and  density  D, 

F  is  the  product  of  the  volume,  density,  and  acceleration  due  to  gravity, 

£  3 

i.e.,  F  =  ( 4 / 3 ttR  )pg.  The  drag  force  F  for  viscous-dominated 
S  v 

Stokes  flow  is  F  =  bTTpRV  where  U  is  the  viscosity  of  the  fluid  and 

Y  oo 

V  is  the  velocity  of  the  particle  through  the  fluid. 

Equating  these  two  forces  and  solving  for  the  velocity  V  ,  we  obtain 


V 

OO 


=  W2 

3  U 


(B-l) 


For  particles  of  given  density  in  a  fluid  of  known  viscosity,  VM  depends 

2 

onLy  upon  R*" .  This  relationship  is  shown  by  the  heavy  solid  line  in  Fig. 


Relationship  of  Particle  Size  to  Settling  Velocity 
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PARTICLE  SIZES  PREDICTED  FROM  COLLISIONLESS  GAS  DYNAMICS 


For  this  case,  the  drag  coefficient  C__  =  F_/l/2p  V4R 
22  D  D  00  00 

is  given  by 


2S2  ♦  I  -S2  4S4  +  4S2  -  l  e/#1,  2(1-0  Ttl/2/Tw 

CD  ■  771^  •  *  - *«»>  *  — h — [~m 


where  the  speed  ratio  S  is  given  in  terms  of  and  the  gas  constant  R  and 
temperature  T  by  S  =  V  //2RT  .  In  the  last  term,  T  is  the 

1  CO  oo  oo  W 

parcicle  temperature  and  £  is  tne  fraction  of  molecular  collisions  with 
particles  that  are  specular  rather  than  diffuse.  In  the  following,  we  assume 
T  =  T  and  £  =  0. 

w  oo 

For  the  settling  velocity  case  in  which  we  are  interested,  <<  2RToo> 
i.e.,  S  is  small.  When  we  substitute  asymptotic  expressions,  valid  for  small 
-S2 

S,  in  place  of  e  and  the  error  function  erf(S),  we  may  combine  the  above 
expressions  to  obtain* 


Equating  this  force  with  the  gravity 
V  we  obtain 

OO 


3 

force  F  =  (4/3  7iR  )pg  and  solving  for 
g 


Comparing  this  with  F,q.  (B-l),  we  see  that  for  the  collisionless  case 
the  dependence  on  R  is  linear  and  depends  inversely  on  the  ambient 
pressure  p  .  For  a  room  temperature  ambient,  Eq.  (B-2)  becomes 
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( B-3  ) 


Vm  =  0.0591  — 

oo  n 

r  oo 


for  in  ft/sec,  R  in  Mm,  and  in  torr.  Representative  results 
calculated  using  Eq.  (B-3)  are  shown  by  the  lighter  lines  in  Fig.  (B-l). 

Thus,  we  conclude  that  for  the  range  of  settling  velocities  seen  in  our 
experiments,  i.e.,  0.011  to  0.033  ft/sec,  the  predicted  particle  sizes  are 

2-4  uni  (for  continuum  flow) 

-2 

8-20  iwi  (for  collisionless  flow,  p0  *  5  x  10  torr) 

-2 

For  =  5  x  )0  torr,  the  mean  free  path  for  molecules  is  approximately 
1000  Mm,  i.e.,  1  mm.  Since  this  is  much  largar  than  the  predicted  continuum 
particle  sizes  of  2-4  p»ii,  the  calculated  particle  sizes  (continuum)  are  not 
likely  to  be  correct,  because  the  flow  is  not  continuum.  Conversely,  the 
particle  sizes  predicted  by  the  collisionless  equations  may  be  somewhat  too 
smal  l . 
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